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Abstract. Mcsoscopic hydrodynamic equations arc solved to investigate the 
dynamics of nanodroplets positioned near a topographic step of the supporting 
substrate. Our results show that the dynamics depends on the characteristic length 
scales of the system given by the height of the step and the size of the nanodroplets 
as well as on the constituting substances of both the nanodroplets and the substrate. 
The lateral motion of nanodroplets far from the step can be described well in terms 
of a power law of the distance from the step. In general the direction of the motion 
depends on the details of the effective laterally varying intermolecular forces. But for 
nanodroplets positioned far from the step it is solely given by the sign of the Hamakcr 
constant of the system. Moreover, our study reveals that the steps always act as a 
barrier for transporting liquid droplets from one side of the step to the other. 
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1. Introduction 

Understanding the wetting behavior of liquids on solid substrates [1, 2] is a prerequisite 
for making use of a myriad of biological and technological applications such as eye 
irrigation, cell adhesion, tertiary oil recovery, coating, lubrication, paper industry, 
micro-mechanical devices, and the production of integrated circuits. Generically, the 
solid surfaces in the above mentioned examples are not ideal in the sense that they 
are neither smooth nor homogeneous. Most surfaces are topographically or chemically 
heterogeneous. Such heterogeneities may substantially change the wetting behavior 
of these surfaces [3], which is not necessarily detrimental with respect to envisaged 
applications. Certain topographically structured surfaces are superhydrophobic or 
superhydrophilic. In the first case droplets roll off these substrates (instead of flowing), 
such that these surfaces are self-cleaning [4, 5, 6, 7, 8, 9, 10, 11]. In the second case 
the surface topography leads to a complete spreading of droplets [12, 13, 14]. Tailored 
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topographic surface structures can induce particular dewetting processes which in turn 
can be exploited to pattern substrates on the micron scale [15, 16]. 

Microfluidics is another strong driving force for the research on the dynamics of 
fluids on structured substrates. Shrinking standard laboratory setups to a lab-on-a-chip 
promises huge cost reduction and speed-up [17, 18]. Open microfluidic systems, i.e., 
with free liquid-vapor or liquid-liquid interfaces, may provide various advantages such 
as reduced friction, better accessibility of the reactants, and reduced risk of clogging 
by solute particles [19, 20, 21, 22]. In open microfluidic devices fluids are guided along 
chemical channels [3, 23, 24] or in grooves [25], which can be chemically patterned in 
oder to provide additional functionality [22]. 

Wetting phenomena on topographically structured substrates have attracted 
substantial research efforts [25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35] with, however, 
the main focus on equilibrium phenomena. In view of the aforementioned applications, 
dynamical aspects are of particular interest. In spite of this demand, theoretical work 
on the dynamics of liquid films and droplets on topographically structured substrates 
has started only recently. In most of these studies the dynamics of the fluids is assumed 
to be well described by macroscopic hydrodynamic equations, which are solved either 
directly [36], by a lattice Boltzmann method [37, 38], or in the thin film (lubrication) 
regime [39, 40, 41, 42]. The applicability of this latter method is limited because 
the inherent long-wavelength approximation does not keep track of many relevant 
microscopic features [43]. 

On the nanoscale, macroscopic hydrodynamic equations turn out to be inadequate 
for describing the dynamics of fluids. Overcoming this deficit is the focus of a 
new research area called nanofluidics [44, 45]. Wetting phenomena in particular 
reveal these deviations; for a recent review of these issues see Ref. [46]. However, 
hydrodynamic equations can be augmented to include hydrodynamic slip, the finite 
range of intermolecular interactions, and thermal fluctuations. The resulting mesoscopic 
hydrodynamic equations have been rather successful in analyzing, e.g., the dynamics 
of dewetting on homogeneous substrates [47, 48]. The presence of intermolecular 
interactions can be summarized into the so-called disjoining pressure (DJP), II = 
— <9$ / dy where the effective interface potential $ is the cost of free energy to maintain 
a homogeneous wetting film of prescribed thickness y. On a homogeneous substrate 
$ is independent of lateral coordinates parallel to the substrate surface and the 
equilibrium wetting film thickness yo minimizes $(?/). However, on chemically or 
topographically inhomogeneous substrates (structured, rough, or dirty) the generalized 
disjoining pressure does depend in addition on these lateral coordinates. In most studies, 
the lateral variations of the disjoining pressure have been modelled rather crudely, 
i.e., the substrate is assumed to be locally homogeneous and lateral interferences of 
heterogeneities are neglected: e.g., a step is typically modelled by an abrupt change of 
the disjoining pressure [49, 50, 51, 52, 53]. 

Recently we have demonstrated, that the actually smooth variation of the lateral 
action of surface heterogeneities can change the behavior of droplets in the vicinity 
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of chemical steps [54, 55] or topographical features (edges and wedges) [56] even 
qualitatively. In the present study we extend these results to the case of an isolated 
straight topographic step in an otherwise homogeneous substrate (as shown in Fig. 1) 
and we recover the previously studied case of isolated wedges and edges in the limit 
of infinite step height h. We should emphasize that our investigation provides only a 
first but nonetheless essential step towards understanding the dynamics of droplets on 
arbitrarily structured substrates. Although more refined than previously used models 
the present one is still rather simple. We only consider additive Lennard- Jones type 
intermolecular interactions, i.e., we do not take into account electrostatic interactions 
which would be very important for polar fluids. We assume the fluid to be Newtonian, 
non- volatile, and incompressible (which is compatible with the frequently used so-called 
sharp- kink approximation of classical equilibrium density functional theory (see, e.g., 
Ref. [57]). We also assume a no-slip boundary condition at the solid surface [58] and 
neglect the influence of thermal fluctuations [59]. For numerical reasons we restrict our 
investigation to two-dimensional (2D) droplets, corresponding to three-dimensional (3D) 
liquid ridges (or rivulets) which are translationally invariant in the direction parallel to 
the step; nonetheless we expect our results to hold qualitatively also for 3D droplets. 

2. Summary 

We study the dynamics of non- volatile and Newtonian nanodroplets (corresponding to 
three-dimensional ridges which are translationally invariant in one lateral direction) on 
topographically stepped surfaces within the framework of mesoscopic hydrodynamics, 
i.e., by solving the augmented Stokes equation presented in Sec. 3 with the numerical 
method described in Appendix A. We consider in particular the effects due to the long 
range of Lennard- Jones type intermolecular interactions which enter the theoretical 
description in terms of the disjoining pressure (DJP) as illustrated in Fig. 1. We assume 
the substrate to be chemically homogeneous in the lateral directions and the surface to 
be covered by a thin layer of a different material. As detailded in Sec. 4 this leads to two 
adjustable parameters B and C which enter into the DJP and characterize the wetting 
properties of the substrate, i.e., the equilibrium contact angle 9 eq and the wetting film 
thickness yo (see Fig. 2). As shown in Fig. 3 both for positive and for negative Hamaker 
constants one can find a one-parameter family of pairs (B, C) leading to the same 9 eq 
on a flat substrate (i.e., without a step). As shown in Fig. 4 nanodroplets on substrates 
with the same 9 eq but with different values of B and C assume shapes which differ 
mainly in the vicinity of the three-phase contact line with the apex region is almost 
unaffected by the substrate potential. 

The results of the numerical solution of the mesoscopic hydrodynamic equations are 
presented in Sec. 5. In contrast to macroscopic expectations based on a capillary model 
(i.e., taking into account only interface energies and neglecting the long range of the 
intermolecular interactions), topographic steps do influence droplets in their vicinity: on 
substrates with a positive Hamaker constant (Figs. 9 and 15), droplets move in uphill 
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Figure 1. Nanodroplets positioned near a topographic step of height h (on the top 
side and bottom side of the step) are exposed to the vertically and laterally varying 
disjoining pressure II, the contour plot of which is shown. The topographic step and 
the drops are taken to be translationally invariant along the z axis (i.e., orthogonal 
to the image plane). In (a) the substrate is chosen to correspond to the minus case 
~] with (B = 0, C = 1) and in (b) the substrate corresponds to the plus case | + | 
(B = —2.5, C = 1) (see Eq. (22) for definitions). Lengths (x, y, h) and the disjoining 
pressure II are measured in units of b and 7/6, respectively (see the main text for 
definitions). 
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direction while on substrates with a negative Hamaker constant (Figs. 5 and 11) the 
droplets move in the opposite direction. As expected the forces on the droplets and 
their resulting velocity increase with the step height, but also with the absolute value of 
the Hamaker constant. This is the case if the contact angle is varied (see, e.g., Figs. 7 
and 12) and even if the contact angle is fixed by varying the Hamaker constant and the 
properties of the coating layer together (see Figs. 8, 10, 13, and 16). 

The speed of the droplets increases with their size as demonstrated in Fig. 14. As 
detailed in Subsec. 6.2, the influence of the step on a droplet can be phrased in terms 
of an effective wettability gradient, i.e., a spatially varying equilibrium contact angle. 
The driving force on droplets on such substrates increases linearly with the droplet size 
because the difference in equilibrium contact angle at the two contact lines of the liquid 
ridges increases roughly linearly with the distance from the steps. 

The velocity of droplets driven away from the step decreases rapidly with the 
distance from the step as shown in Sec. 6. But droplets moving towards the step (either 
on the top side or on the bottom side of the step) stop with their leading contact line 
close to the step edge or wedge, respectively. Therefore they do not cross the step (see 
Figs. 5, 6, 7, 8, 15, and 16). Accordingly, edges, wedges, and steps act as barriers for 
migrating droplets (which is also true macroscopically) because droplets sitting right 
at the tip of an edge are in an free-energetically unfavorable state (see Fig. 17) while 
droplets located in the corner of a wedge are in a state corresponding to a local minimum 
of the free energy (see Fig. 20). Therefore, an external force is required to push droplets 
over edges (see Fig. 18) or to pull them out of wedges (see Fig. 20). In both cases, 
the total (i.e., integrated over the droplet volume) force required to accomplish this 
increases slightly with the droplet volume, but less than linearly. This means, that 
if the force is applied via a body force density acting per unit volume (e.g., gravity) 
larger droplets experience a larger force and therefore overcome steps more easily. In 
addition, the lateral action of intermolecular forces can also pin droplets at edges and 
near wedges. However, droplets which initially span a topographic step always end up 
filling the wedge at the step base, either with the upper contact line pinned at the step 
edge or, if the droplet volume is too small, with the upper contact line on the vertical 
part of the step, as shown in Fig. 19. 

A deeper understanding of the dynamics of droplets in the vicinity of edges and 
wedges can be reached by analyzing the forces acting on the droplet surface, i.e., 
the disjoining pressure and surface tension (see Eqs. (6) and (7), respectively). As 
demonstrated in Fig. 6, if the droplets move under the influence of the topographic step 
only, the main contribution to the driving force stems from the disjoining pressure. As 
shown in Figs. 21 and 23, the numerically observed features of the dynamics of droplets 
can be understood in terms of the disjoining pressure induced force density on the 
droplets calculated for droplets of simple parabolic shapes used as initial conditions for 
the numerical solution of the hydrodynamic equations. As shown in Fig. 22 the actual 
relaxed droplet shape is different but the calculated forces depend only weakly on the 
deviation of the actual shape from its parabolic approximation. In the limit of large 
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distances from the step the force can be calculated analytically (see Subsec. 6.2): far 
from the step the total force per unit ridge length F n = fn A d (with the cross-sectional 
area Ad) essentially depends on the ratio of the step height h and the distance from the 
step x as well as on the ratio of the apex height y m and x. The corresponding asymptotic 
results are summarized in Fig. 24. In all cases the force density varies according to a 
power law x~^ with £ e {3,4,5}. For finite sized droplet and steps of finite height we 
obtain the fastest decay and for almost macroscopic droplets in the vicinity of finite sized 
steps as well as for nanodroplets near isolated edges and wedges we get ( = 4. While 
our present analysis cannot be applied to the case of an almost macroscopic droplet 
in a wedge, for large drops (y m /x — > oo) next to an isolated edge we get the weakest 
decay with ( = 3. In any case, the total force per unit length Fn is proportional to 
the Hamaker constant as observed in the numerical solution of the mesoscopic Stokes 
dynamics as well as in the force analysis presented in Subec. 6.1. The dynamics of large 
drops (y m /x — > oo) is equivalent to the dynamics of macroscopic drops on a surface with 
an effective chemical wettability gradient (i.e., a spatially varying "equilibrium contact 
angle" 6 eq (x)) [60, 61, 62]. 



3. Mesoscopic hydrodynamic equations 

At low Reynolds numbers the mean field dynamics of an incompressible Newtonian fluid 
of viscosity /i is given by the Navier-Stokes equation for the local pressure p(r, t) and 
the flow field u(r,t): 

V • <x = -Vp + fi V 2 u = 0, (1) 
V-u=0, (2) 

with the stress tensor = —pSy + /i (c^-Uj + diUj). In this study, we neglect the 
influence of the vapor phase or air above the film. Therefore the tangential components 
of the component of the stress tensor cr • n normal to the liquid- vapor surface Ti v (with 
outward pointing normal vector n) is zero. The normal component of cr ■ n, i.e., the 
normal forces acting on the liquid surface, are given by the sum of the Laplace pressure 
and of the disjoining pressure: 

cr-n = n( r yK + U + xg) at Ti v , (3) 

with the surface tension coefficient 7 and the local mean curvature k of the liquid 
surface; g is the strength of a spatially constant external body force density pointing 
in the x-direction (with — gx as the corresponding potential) which we introduce in 
order to study the strength of barriers to the lateral motion of droplets. Alternatively, 
for incompressible fluids one can define a new pressure p' = p — xg such that the 
external body force density g enters into the Stokes equation Eq. (1) rather than the 
boundary condition in Eq. (3): = — Vp' + e^ g + yuV 2 u . Although this approach might 
be more intuitive, the equivalent form used here is more convenient for implementing 
the boundary element method used here to numerically solve these equations (see 
Appendix A). 
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The dynamics of the free liquid surface is determined by mass conservation together 
with the incompressibility condition: the local normal velocity is identical to the normal 
component of the local flow field. 

We neglect hydrodynamic slip at the liquid-substrate surface Ti s and we only 
consider impermeable substrates. Since we assume the substrate to be stationary this 
results in the following boundary condition for the flow field: 

u = at T ls . (4) 

In order to avoid strong initial shape relaxation of the droplets (in response to 
placing them on the substrate with a certain shape) which can lead to significant lateral 
displacements [63], we choose a parabolic initial profile which is smoothly connected to 
a precursor film of thickness yo'- 

/I -IN 2' 

/ IX — XI \ 

1 - 



y(x\ t = 0) = y + a 



(5) 



such that a is the droplet height at the center and half the base width. Accordingly the 
distance of the droplet edge from the step at x = is given by £ = \x\ — a with x the 
position of the center of the droplet in the x-direction. The parameter m specifies the 
smoothness of the transition region from the drop to the wetting layer. In this study 
we choose m to be 10. We investigate the droplet dynamics for two different situations. 
In the first one we position the droplet on the top side of the step of height h with the 
three-phase contact line (x = x + a, y = h + yo, z) at a distance I = —x — a with x < —a 
from the step edge at x — 0. In the second situation we place the droplet on the bottom 
side of the step with the three-phase contact line (x = x — a, y = y , z) at a distance 
t = x — a with x > a from the wedge at the base of the step. 

In equilibrium Eq. (3) reduces to the Euler-Lagrange equation of the effective 
interface Hamiltonian of a fluid film on a substrate as derived, e.g., in Ref. [64]. 
This means that we approximate the normal forces on the liquid surface due to the 
intermolecular interactions by the disjoining pressure derived for equilibrium systems. 

In a non-equilibrium situation, the unbalanced forces acting on the fluid surface 
add up to a resulting net force on the liquid body. We separately consider the two 
contributions /n and f 1 from the disjoining pressure and from the Laplace pressure, 
respectively, both normalized by the droplet volume Qd and given by the following 
integrals over the liquid-vapor surface of the droplets: 

/n(x) = — / n(x, y) n x ds (6) 



d 

fi( x ) = FT / lKn x ds. (7) 

For a liquid ridge translationally invariant in ^-direction both integrals as well as 
Vtd = AdL are proportional to the macroscopic ridge length L, so that the latter drops 
out of the expressions for the force densities (in units of N/m 3 ) fu and / 7 . In three 
dimensions ds is a two-dimensional surface area element. Ad is the two-dimensional 
cross-sectional area of the liquid ridge. 
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4. Model of the heterogeneity 

In the following we calculate the disjoining pressure for a fluid film or droplet near 
a topographic step as displayed in Fig. 1. Apart from a very thin coating layer of 
thickness d we assume the substrate material to be homogeneous, disregarding its 
discrete molecular structure. Many substrates used in experiments are coated, e.g., 
by a native oxide layer or by a polymer brush which is used to modify the wetting 
properties of the substrate. However, a more refined analysis of the DJP, which takes 
the molecular structure of the substrate and of the fluid into account, yields terms of a 
form similar to those generated by a coating layer [2, 65]. In general, i.e., far from the 
critical point of the fluid, the vapor or gas phase covering the system has a negligible 
density which we neglect completely. Assuming pairwise additivity of the intermolecular 
interactions, i.e., the fluid particles as well as the fluid and the substrate particles are 
taken to interact with each other via pair potentials V a p(r) where a and (3 relate to 
liquid (I), substrate (s), or coating (c) particles and r is the interatomic distance, one 
can show that the disjoining pressure (DJP) of the system is given by [26] 

n(r) = f [p 2 V u (r - r') - Pl p s V sl (r - r')] d\ (8) 
Jn s 

with r, r' £ and p\ and p s as the number densities of the liquid and substrate, 
respectively. Q s is the actual substrate volume. 

In order to facilitate the calculation of the disjoining pressure of the step we 
decompose it into contributions from quarter spaces (edges) forming building blocks 
which can be calculated analytically. We first consider an edge occupying the lower left 
quarter space ttf = {r£]R 3 |x<0Ai/<0}, which in the following we denote by <. 
For Lennard- Jones type pair potentials V a p(r) = M a p/r 12 — N a p/r 6 , where M a p and 
N a p are material parameters, the DJP in the vicinity of a non-coated edge occupying 
flf is given by 

-oo 

n e «(*, y)= j dx' J dy> J dz> (j^p " r^je) , (9) 

— oo — oo — oo 

where AM e = pf Mu — pi p s M\ s and AN e = pf Nu — pi p s Ni s . The first term dominates 
close to the surface of the edge and the second term at large distances from the substrate. 

All integrals in Eq. (9) can be calculated analytically and one obtains the DJP as 
the corresponding difference = AM e J e 12<1 — AN e 1^ of two contributions with 

W*>4 = 11520 xVW) 7 / 2 h28 ° XW + 
y 4 ) - U8x 2 y 2 (x 12 + y 12 ) - 128 (x 16 + 

y 16 ) + 128 (x 9 + y 9 )(x 2 + y 2 ) 7 / 2 

- 35 x 8 y s - 560 x 4 y 4 (x 8 + y 8 )] 

(10) 
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and 

y) = r— r-p 7====[ 2 (* 3 + y 3 ) V^TV 2 

-2{x A + y 4 )-y 2 x 2 }- (11) 

The contributions to the disjoining pressure of a thin coating layer of thickness d 
on the -upper side of the edge occupying = {r 6 I 3 | x < 0, —d < y < 0}, the right 
part of the edge occupying = {r 6 R 3 —d < x < 0, y < 0}, and the thin rod 
which fills the tip area of the edge = {r G I 3 — d < x < 0, —d < y < 0} can be 
calculated analogously: 

Jo*** |r - r I |r - r| 

(12) 

with A M c = pfMu - p c pi M d and A N c = pf N u - p c p x N d ; \ stands for u (wpper), r 
(right), or t (tip). Actual coating layers have a more complicated structure, in particular 
in the direct vicinity of edges and wedges, which depends on the specific combination 
of coating and substrate material as well as on the way the coating is produced. Such 
details can influence droplets if their contact line is right at the edge or wedge but the 
effect is proportional to the square of the coating layer thickness d. For simplicity we 
only consider systems with coating layers which are thin compared to the wetting film 
thickness (see below), for which the contribution from the thin rod of coating material at 
the tip of the edge or in the corner of the wedge is irrelevant. According to Eq. (12) the 
contribution to the disjoining pressure from the upper coating layer can be decomposed 
into II" = AM/ C 12n (x, y) — AN I® u (x,y). To first order in d we obtain 

1280 (x 2 + y 2 ) 9/2 y 10 

[128 (x 2 + y 2 ) 9/2 - 315 x y 8 - 840 x 3 y 6 - 

1008 a; V - 576 x V - 128 x 9 ] (13) 

Ti d 



and 



,4 (rnl _L „,2^ 3 / 2 



\y* {x z + y z 

[-2(x 2 + y 2 f 12 + 3xy 2 + 2x 3 ]. (14) 

By symmetry one has IF^fa, y) = H'^ < (y,x) for the contribution of the vertical part of 
the coating. The DJP of a coated edge occupying fife = {r G IR 3 | x < A y < 0} is 
therefore given by 

n«(2, y) = Uf(x + d,y + d) + I%<(x, y) + U r <(x, y). (15) 

The DJP contribution from a coated edge occupying the right quarter space 
^ce = {r 6 R 3 | ^ ) A i/ < 0} can be obtained analogously. However, since the 
integrals for the right part corresponding to Eqs. (9) and (12) are the mirror image 
(with respect to the yz-plane) of their counterparts for the left hand side, the former 
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ones can be expressed in terms of the latter ones. Therefore the D JP of the coated lower 
right quarter space n^(x, y) is equal to n^,(— x, y). Combining the contributions of the 
left and the right part leads to the following expression for the DJP of a step of height 
h: 

U(x, y) = U<{x, y + h) + U<(-x, y) - 2 IT<(x, y). (16) 

The last term on the right hand side of Eq. (16) removes the artificial extra coatings on 
the left and the right quarter spaces (at x = 0, y < 0) which get buried upon building 
the step out of the coated edges. Figure 1 shows typical examples for the DJP. The DJP 
is not only a function of the vertical distance from the substrate, but also of the lateral 
distance from the step. In this regard, the substrate in the vicinity of the step resembles 
a chemically structured substrate with laterally varying wettability [60, 61, 62]. 

For positions far from the step the distribution of the DJP resembles that of the 
coated, laterally /iomogeneous flat substrate obtained by setting h = in Eq. (16). To 
linear order in d one has 

vr AM e ir AN e vr AM e d n AN e d 
<*W) = ~J^y~9 ^3 5^10 + 2y 4 

nAM c d nAN c d 
5y 10 2y* ' 1 j 

Since the repulsive contributions decay rapidly with distance from the substrate we 
neglect all those repulsive contributions which are shorter ranged than the corresponding 
term (~ y~ 9 ) arising from Ilg 2<I (x, y) [57, 65], leading to 
^ nAM e 7rAiV e n AN d 

45 y y Qy A 2 y 4 
with AiV = AN C — AN e . The equilibrium thicknesses yo of the wetting film on such a 
substrate minimizes the effective interface potential [65, 66] 

POO 

$ ch (y) = / Ii ch {y)dy. (19) 



With Eq. (18) this leads to 

. , vrAM e vrAiV e -n AN d , , 

®ch(y) = -^rr ~ ^r~2 ^-3-- 20 

360 y 6 12 y z by 6 
The second term is usually written as —H e / (12 ny 2 ), where H e = tt 2 AN e is the so-called 
Hamaker constant. 

At this point we introduce dimensionless quantities (marked by *) such that lengths 
are measured in units of b = [2|AM e |/(15 |AiV e |)] 1/6 which for AM e > and AN e > 
is the equilibrium wetting film thickness y^ on the uncoated flat substrate. The DJP is 
measured in units of the ratio 7/6 where 7 is the liquid- vapor surface tension. Thus the 
dimensionless DJP U*j- = U c f 6/7 far from the edge has the form 
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In the first and second term of Eq. (21) the upper (lower) sign corresponds to AM e < 
(AM e > 0) and AiV e < (AN e > 0), respectively. The dimensionless amplitude 
C = Ab/j, with A = 7r(|AM e |/45)~ 1/2 (|AiV e |/6) 3/2 , compares the strength of the 
effective intermolecular forces in the uncoated case and of the surface tension forces. 
The amplitude B = 7iANd/(2Ab A ) measures the strength of the coating layer. Since 
the molecular structure of the substrate and of the fluid yields a term of the same form 
[2, 65] we consider B itself as a parameter independent of the actual properties of the 
coating layer. For the interactions considered here, AM e > is a necessary condition 
for the occurrence of an equilibrium wetting layer of nonzero thickness but AN e can be 
positive or negative. Therefore the first term in Eq. (21) can only be positive while the 
second term can be positive or negative. In the following we shall refer to these two 
cases simply as the minus ( | — | ) and the plus ( | + 1 ) case. In order to avoid a clumsy 
notation in the following we also drop the stars. With this, one has 

n*(y) = c ( - 9 t - z + 4) • (22) 
\y 9 y y J 

Figure 2 shows the typical profile of U c h(y) f° r the minus and the plus case and also 
the corresponding equilibrium wetting layer thickness yo for which U cn (yo) = 0. While 
the parameter C measures the strength of the DJP, by changing B one can modify the 
shape of the DJP [56]. In Eq. (22) the admissible value ranges of C and B which provide 
partial wetting can be inferred from considering the equilibrium contact angle 9 [2]: 

poo 

cos 6= 1+ / U ch (y)dy. (23) 

The admissible value ranges of B and C for which 0° < 9 < 180° (partial or incomplete 
wetting) are given in Fig. 3 for both the minus and the plus case. In the minus case, for 
each value of B one can find a value of C such that the resulting substrate is partially 
wet. Since the signs of the first two terms in Eq. (22) differ the disjoining pressure has 
a zero for any B and the depth of the minimum of the corresponding effective interface 
potential can be tuned by choosing an appropriate value for C. In the plus case, however, 
B has to be negative in order to obtain a sign change of II. The maximum admissible 
value of B (i.e., B < B max ) can be obtained by simultaneously solving the following 
equations for y and B max : 

n ch (y ) = ^-\ + ^P = (24) 
Vo Vo Vo 

**(ifo) = A- A + irr = o, (25) 

%o 2 Vo %o 
from which one finds B max = —1.868 (compare Fig. 2(d)). 

In order to obtain dimensionless hydrodynamic equations (see Eqs. (l)-(3)) we 
choose Ab/fi as the velocity scale. With this, the dimensionless form of the stess tensor 
is given by = —pSij + C (djUi + diUj) and the surface tension coefficient drops out 
of Eq. (3). The dimensionless time is given in units of fi/A. In order to study the 
dynamics of nanodroplets we solve the dimensionless hydrodynamic equations with a 
standard biharmonic boundary integral method described in more detail in Appendix A. 
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Figure 2. Typical (B = —2) DJP (in units of C) of a fiat homogenous substrate for 
(a) the minus and (b) the plus case (see Eq. (22)). The corresponding zeros yo of the 
DJP for different values of B are given in (c) and (d) for the minus and the plus case, 
respectively. In (c) and (d) full lines indicate stable wetting films and dashed lines 
unstable films. 




Figure 3. The value ranges of B and C for which the system exhibits a partial 
wetting (PW) situation, i.e., 0° < 9 eq < 180° for the minus (a) and the plus (b) case. 
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Figure 4. Top panel: Droplet prepared in its initial configuration. Bottom panel: 
The equilibrium configuration of a nanodroplet with a = 15 on a flat homogeneous 
substrate for various values of B and C for the minus and the plus case. (B, C) for 
the minus case are I (-1, 7.7583) with y = 0.88, II (0, 2.6667) with y = 1, and III 
(1, 1.2703) with y = 1.3, and for the plus case I (-2.5, 4.2327) with y Q = 0.91 and II 
(-4, 0.9265) with j/q = 0.79. The values of B and C are chosen such that in all cases 

e eq = 90°. 



5. Results 

5.1. N anodroplets on homogeneous flat substrates 

In order to provide the information and terminology required for the subsequent 
considerations we first recall some basic results for the wetting of flat and homogeneous 
substrates. For this purpose a nanodroplet with a = 15 and an initial configuration 
given by Eq. (5) was positioned on the substrate. Figure 4 shows the equilibrium profile 
of the nanodroplet for various values of C and B resulting in an equilibrium contact 
angle 9 eq = 90° for both the minus and the plus case, i.e., the values (B,C) lie on the 
dashed curves in Figs. 3(a) and 3(b). It is evident from the figure that the droplets have 
relaxed from the initial condition. The equilibrium profiles in all cases are roughly equal 
but the nanodroplets differ near their contact lines (see the inset of Fig. 4) and with 
respect to their heights. The term proportional to B in Eq. (22) is rather short-ranged 
and most important in the direct vicinity of the substrate. The top parts of the droplets 
are only influenced by the term C/y 3 such that the curvature at the peak changes with 
C, independently of B. This also changes the droplet height. However, also the wetting 
film thickness y® changes with B, such that the differences in droplet height in Fig. 4 
are a combined result of both effects. 

Due to the translational symmetry of the substrate and due to the symmetry of the 
initial drop configuration the shape relaxation does not result in a lateral displacement 
of the droplets, in contrast droplets placed on heterogeneous substrates [63, 67]. 
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Figure 5. (a) The effect of the step height on the dynamics of nanodroplets positioned 
on the top side of a step for the minus case. The dashed and the solid lines correspond 
to nanodroplets at t = 200 and t = 18700, respectively. The droplets arc initially 
positioned at a distance I — 10 from the step. The droplets have an initial radius 
a = 15. C = 2.667 and B = correspond to 8 eq = 90°. The vertical scale is equal to 
the lateral scale. The corresponding lateral (b) and vertical (c) position of the center 
of mass (cc, y) of the droplet relative to the step edge as a function of time. 

5.2. Nanodroplets on the top side of steps 

Previous studies of droplets near edges (corresponding to steps of infinite height) have 
shown that, in contrast to what is expected from a simple macroscopic model taking into 
account only interface energies, droplets are attracted towards the edge in the minus 
case and repelled from the edge in the plus case [56]. In the minus case, the droplets 
move towards the edge with increasing velocity, but they stop rather abruptly before the 
leading contact line reaches the edge. The distance from the edge at which the droplets 
stop increases with decreasing B, i.e., with increasing strength of the coating layer. In 
the plus case, the droplets move away from the step with a velocity which decreases 
with the distance from the step. The strength of the attraction or repulsion is expected 
to be lower for steps of finite height. 

5.2.1. Minus case In order to test the influence of the step height on the dynamics 
of nanodroplets identical droplets of half base width a = 15 were placed at a distance 
i = 10 from steps of height h =2.5, 5, 10, 15, 20, and oo. The results of the numerical 
solution of the mesoscopic hydrodynamic equations for the minus case are shown in 




U i i i I i i i I i i i L 
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x s 

Figure 6. (a) The initial and the final configurations of a droplet with a = 15. Initially 
it is positioned at £ = 10 (x = —25, dash-dotted line). For the minus case it moves 
towards an isolated edge where it stops with the leading contact line pinned at the 
step edge (full line). For this final configuration the arclength s of the interface is 
measured as indicated from a certain position on the wetting layer on the vertical side 
of the step, (b) A comparison between DJP induced (Eq. (6), full line) and surface 
tension induced (Eq. (7), dashed line) lateral force densities during the motion (with 
the leading three-phase contact line still well separated from the edge) expressed in 
terms of the position x of the center of mass of the droplet, (c) Laplace pressure 
7 k (dash-dotted line) and DJP II (full line) on the surface of the droplet in the final 
equilibrium configuration as a function of the arclength s. In the absence of other 
external forces (e.g., g in Eq. (3)) at each point on the droplet surface these add up 
to the constant pressure p inside the droplet. B = and C — 2.667 correspond to 
9 eq = 90°. Force densities (b) and pressures (c) are measured in units of 7/b 2 and 7/6, 
respectively. 
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Figure 7. The effect of the contact angle on the dynamics of droplets on the top 
side of the step for a droplet with initial height a = 15 and for the minus case, (a) 
The initial profile (£ = 10) is shown in the top panel. The lower graphs show the 
configurations of the droplets after the initial relaxation (t/C « 60, dashed lines) and 
in the final stages (t/C =6200, 7400, and 7800 for I, II, and III, respectively, solid lines) 
for 9 eq = 75.5° (I), 97.2° (II), and 120° (III) from top to bottom. The corresponding 
substrate parameters are I (C = 2, B = 0), II (C = 3, B = 0), and III (C = 4,B = 0), 
respectively. In (b) and (c) as function of time the corresponding lateral and vertical 
positions x and y, respectively, of the center of mass of the droplets are shown relative 
to the step edge. The dips in (c) occur when the leading three-phase contact line 
reaches the edge; then the droplet stops (compare to (b)). 



Fig. 5(a) for C = 2.667 and B = which corresponds to 9 eq = 90°. In order to have a 
better view on the dynamics we monitor the time evolution of the position of the center 
of mass of the droplets (x, y) relative to the step edge in Figs. 5(b) and 5(c), where x 
and y are given by 

L x dV f n y dV 

* J^v' S = J^v- h ' (26) 

with Qd denoting the droplet volume. Since the droplets are smoothly connected to the 
wetting film, which on large substrates would influence the center of mass of the fluid, 
in calculating x and y we only consider the fluid above y = Coyo + h with Co > 1, i.e., 
only the fluid volume slightly above the wetting film. We selected c = 1.2; but since 
we focus on substrates with equilibrium contact angles of about 90° the results are only 
weakly affected by the precise choice of the value of c . 
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Figure 8. The effect of changing B (while varying C such that 9 eq = 90° in all cases, 
see Fig. 3(a)) on the dynamics of the droplets on the top side of an edge for an initial 
droplet with a = 15 and for the minus case. The values of B and C arc I (B = — 1, 
C = 1.2703), II (B = 0, C = 2.6667), and III (B = 1, C = 7.7583). (a) The top 
panel depicts the initial profile (£ = 10). The dashed lines show the configuration of 
the droplet after the initial relaxation t/C « 60 and the solid lines correspond to the 
final configuration of the droplets at t/C = 16400, 7300, and 2500 for I, II, and III, 
respectively. In (b) and (c) as a function of time the corresponding lateral and vertical 
positions x and y, respectively, of the center of mass of the droplets arc shown relative 
to the step edge. 



In all cases the dynamics of the droplets proceeds in three stages. The first stage is 
a fast initial shape relaxation, similar to the behavior on homogeneous substrates, which 
is accompanied by a lowering of the droplet center of mass y without any considerable 
lateral motion. This is followed by a relatively slow lateral motion towards the edge, 
during which the changes in the droplet shape are almost unnoticeable. Although the 
droplet shape is slightly asymmetric the lateral surface tension induced force density / 7 
defined in Eq. (7) is much smaller than the force density /n induced by the DJP (see 
Eq. (6)) as shown in Figs. 6(b). Figure 5(b) clearly shows that the lateral motion of 
the droplet slows down rapidly as soon as its leading three phase contact line reaches 
the edge. During this third and final stage a part of the droplet volume leaks into the 
wetting film on the vertical part of the step and as a result the droplet experiences a 
sudden drop in its height y (see Fig. 5(c)). The trailing three-phase contact line of the 
droplet still continues its motion towards the step and as a consequence the mean height 
of the droplet increases again and becomes even larger than during the migration stage. 
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While the droplet contracts, its asymmetry gradually increases such that the surface 
tension induced force density / 7 grows and finally, as the equilibrium configuration is 
reached, cancels /n. (This latter stage of cancelation is not visualized in Fig. 6(b) due to 
numerical problems in evaluating the force densities on droplets once they have reached 
the step edge.) In equilibrium, at each point on its surface the Laplace pressure and 
the disjoining pressure ad up to the constant value of the hydrostatic pressure in the 
droplet (see Fig. 6(c)). 

Increasing the step height from h = 2.5 to 5 and 10 results in a significant increase 
in the droplet speed during the migration phase. The asymptotic speed for isolated 
edges (corresponding to h = oo), i.e., the maximum speed, is almost reached for h = 20. 
This height value is large compared to the thickness of the wetting layer but comparable 
with the droplet size; here the base diameter is 2 a = 20. However, in order to be able to 
conclude that the step height above which the droplet perceives the step as an isolated 
edge is comparable with the droplet size further calculations for droplets of different size 
are needed. 

Changing the equilibrium contact angle 6 eq by increasing C while keeping B = 
does not qualitatively change the behavior of the droplets, as shown in Fig. 7(a) for 
droplets with a = 15 near an isolated edge (corresponding to h = oo), apart from the 
increase of y during the initial relaxation process for large 9 eq . The reason for this 
increase in droplet height is, that the initial shape of the droplet is not adapted to the 
substrate parameters. Changing C does not change the functional form of the DJP, 
only its strength. Consequently, droplets move faster for larger C (resulting in larger 
9 eq ) and their final shape is less symmetric. With the leading contact line pinned right 
at the step edge, large 6 eq > 90° also result in an overhang over the step edge. Since 
for fixed wetting film thickness b on the uncoated flat substrate the time scale used to 
obtain dimensionless hydrodynamic equations depends on the substrate parameters in 
the same manner as the dimensionless parameter C, we rescale time by C in Figs. 7(b) 
and 7(c), as well as in all subsequent figures which compare x(t) and y(t) for different 
values of C. This corresponds to changing the substrate material but keeping surface 
tension, viscosity, and wetting film thickness constant [54]. 

Figure 8(a) shows the effect of changing the value of B (while keeping the contact 
angle constant) on the dynamics and on the final configuration of droplets which start 
with a = 15 and for the minus case on the top side of an isolated edge (h = oo). For 
each value of B we choose C such that 9 eq = 90°, i.e., corresponding to the dashed curve 
in Fig. 3(a). For all values of B the droplets move towards the edge. Changing the 
values of B and C does not qualitatively change the behavior of the system. However, a 
closer examination of x and y (see Figs. 8(b) and 8(c), respectively) reveals quantitative 
differences in the dynamics and in the final configuration of the droplets despite the 
fact that the contact angle is the same for all these cases. For larger (positive) values 
of B (and thus C, see Fig. (3) (a)) droplets move faster in lateral direction although the 
contact angle equals 9 eq = 90° for all of them. In addition, the final position of the 
droplets is closer to the step edge for larger values of B, eventually leading to a slight 



CONTENTS 



20 




Figure 9. (a) The effect of the step height on the dynamics of nanodroplets positioned 
on the top side of a step for the plus case. The dashed and the solid lines correspond 
to times t = 200 and t = 10 5 , respectively. The droplets of size a = 15 are initially 
positioned at a distance £ — 10 from the step. C = 4.2327 and B = —2.5 result in 
an equilibrium contact angle 9 eq = 90°. Time evolution of the (b) horizontal position 
x and (c) vertical position y of the center of mass of the droplets relative to the step 
edge. Since y depends only weakly on h only the case h = 20 is shown. 

overhang. The small differences in y for different values of B is related to the fact that 
the shape of nanodroplets is not only determined by eq , as shown in Fig. 4, and that 
the wetting film thickness depends on B. 

5.2.2. Plus case Even if they exhibit the same equal equilibrium contact angles 9 eq the 
behavior of droplets in the plus case differs substantially from that in the minus case: 
the direction of motion is reversed. However, apart from this sign change, the influences 
of the step height, the equilibrium contact angle, and B are similar. 

The dependence of the droplet dynamics on the step height for the plus case is shown 
in Fig. 9(a). The initial size of the droplets is a = 15 and the contact angle 9 eq = 90° 
(with C = 4.2327, B = —2.5). The corresponding lateral and vertical positions of the 
center of mass of the droplets relative to the step edge are shown in Figs. 9(b) and 
9(c), respectively. As in the minus case the migration phase is preceded by a fast initial 
relaxation process (during which y drops slightly). However, the droplets are repelled 
from the step. The lateral speed of the motion continuously decreases as the distance 
of the droplets from the step edge increases. For higher steps the droplets are faster. 
But as in the minus case, the maximum speed (reached for h — oo, i.e., in the case of 
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Figure 10. The effect of B on the dynamics of the droplets near an edge for the 
plus case and for an initial droplet shape with a = 15. The values of B and C, i.e., 
I (B = -2.5, C = 4.2327) and II (B = -4, C = 0.9265) are selected such that the 
equilibrium contact angle is eq = 90° in both cases, (a) The initial distance from the 
step edge is £ — 10. The dashed lines show the droplets after the initial relaxation 
at t/C w 60 and the solid lines correspond to the droplets in the migration phase at 
t/C = 24000. The horizontal position x and the vertical position y of the center of 
mass of the droplet are shown as a function of time in (b) and (c), respectively. For 
less negative values of B the velocity dx/dt is larger (I). 



an isolated edge) is almost reached for h = 20 (see Fig. 9(b)). 

The results for different values of B, while keeping the contact angle 9 eq = 90° 
fixed, are depicted in Fig. 10(a). The corresponding lateral x and vertical y positions of 
the center of mass of the droplet relative to the step edge are given in Figs. 10(b) and 
10(c), respectively. For all the cases considered the droplets move away from the step. 
However, as in the minus case, the droplet speed increases with B (i.e., for less negative 
values of B), even though the contact angle is not changed. The reason for this is, that 
larger (i.e., less negative) values of B require larger values of C in order to maintain the 
same 9 eq . As in the minus case the droplet height, i.e., y, depends on B as well. 

5.3. N anodroplets at the step base 

In Ref. [56] we have demonstrated that droplets near corners, i.e., at the base of a 
step of infinite height, are attracted to the corner in the plus case and repelled from 
the corner in the minus case (while in a macroscopic model taking into account only 
interface energies the free energy of the droplets is independent of their distance from 
the corner). In other words, the direction of motion is reversed as compared to the case 
of the edge. However, as we shall show in the following, at a step composed of an edge 
and a corner at its base, the direction of motion of nanodroplets is the same on both 
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sides of the step. 

5.3.1. Minus case As in the case of droplets on the top side of steps, the step height 
influences the droplet velocity but not the direction of motion and the transition from 
a planar substrate (h = 0) to an isolated wedge (h = oo) is continuous. This is 
demonstrated in Fig. 11 for droplets of size a = 15 starting at a distance i = 10 
from the corner. The initial distance £ is chosen such that after the initial relaxation 
which precedes the migration phase the contact line facing the corner is well separated 
from the wetting layer on the vertical part of the step. For the minus case Fig. 11(a) 
presents the results of our Stokes dynamics calculations for droplets at the step base for 
different step heights. The droplets are repelled from the step and move away with a 
speed which decreases continuously with the distance from the step. The differences in 
droplet speed are significant between step heights h = 2.5, 5, and 10 (see Fig. 11(b)). 
Increasing the step height further influences the dynamics of the droplets only at large 
distances from the step. 

Changing the equilibrium contact angle 6 eq does not change the droplet dynamics 




Figure 11. (a) Nanodroplcts with initial size a = 15 positioned at the base of 
topographic steps of different height for the minus case. The droplets start at a distance 
£=10 from the step. C = 2.667 and B ~ correspond to 6 eq = 90°. The dashed 
and the solid lines correspond to the configurations just after the initial relaxation at 
t = 170 and to a later time t = 30000, respectively. As function of time the horizontal 
position x and the vertical position y of the center of mass relative to the step base 
are shown in (b) and (c), respectively. Since y depends only weakly on h, in (c) only 
the trajectory for h = 20 is shown. 
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Figure 12. (a) Droplets of initial size a = 15 near a corner of substrates with different 
contact angles eq in the minus case. The top panel shows the initial droplet profile with 
£ = 5. The other panels show the droplets after the initial relaxation (t/C = 44.5 (I), 
20.3 (11), and 34.25 (III), dashed lines) and during the migration stage (t/C = 26100 
(I), 17200 (II), and 13600 (III), solid lines) for 6 eq = 75.5° (I, C = 2, B = 0), 97.2° 
(II, C = 3, B = 0), and 120° (III, C = 4, B = 0), respectively. The time evolution of 
the center of mass (x,y) relative to the corner is shown in (b) and (c) for the lateral 
and vertical direction, respectively. 



qualitatively. This is demonstrated in Fig. 12(a) for droplets on substrates with different 
values of C (while keeping B = 0). The top panel shows the initial shape used in all cases 
considered here for the numerical solution of the Stokes dynamics. The corresponding 
lateral position x and vertical position y of the center of mass of the droplets relative 
to the corner are depicted in Figs. 12(b) and 12(c), respectively. Increasing 9 eq (by 
increasing C) results in faster droplet motion. Since the initial droplet shape is not 
adapted to the modified contact angle 6 eq ^ 90°, y changes rapidly during the initial 
relaxation process for 9 eq ^ 90°. 

The dynamics of droplets on substrates with the same contact angle 9 eq = 90° but 
different values of B (with C adapted accordingly) is shown in Fig. 13(a). The top 
panel shows the initial configuration. For all cases considered the droplets move away 
from the step. Comparing x(t) for different values of B (see Fig. 13(b)) shows, that 
the droplet velocity increases with B (which, for fixed 9 eq = 90°, implies increasing C). 
After the initial relaxation process the vertical coordinate y of the center of mass does 
not vary in time (see Fig. 13(c)). 

The droplet dynamics depends on the droplet size. This is demonstrated in Fig. 14 
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for droplets of initial sizes a = 10 and a = 15 starting at x = 20 near an isolated wedge 
for the minus case. The larger droplet moves faster because its two three-phase contact 
lines have a larger lateral distance from each other such that they experience a larger 
difference in the local disjoining pressure. 

5.3.2. Plus case In the plus case the direction of motion of the droplets is reversed as 
compared to the minus case. As shown in Fig. 15(a), the migration speed increases with 
the step height, but the droplets stop before the leading contact line reaches the wedge 
such that the droplets do not move into the corner. As in the other cases discussed so 
far, the droplet speed increases significantly as the step height is increased up to h — 20. 
In Fig. 15(b) the trajectories x(t) for h = 20 and for h = oo almost coincide. The final 
distance of the droplets from the wedge decreases with the step height, but it remains 
finite in the limit h — > oo. Once the droplets reach the wedge there is a brief drop of 
y due to fluid leaking out of the droplet into the corner area. After that their vertical 
position y of the center of mass increases again (see Fig. 15(c)). This increase is the 
result of a contraction of the droplets, which is also observed for droplets on the top 
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Figure 13. (a) Nanodroplets on substrates with different B but the same contact 
angle 9 eq = 90° near a corner in the minus case: (I: B = — 1, C = 1.2703), (II: B = 0, 
C = 2.6667), and (III: B = 1, C = 7.7583), top to bottom. The top panel depicts the 
initial droplet shape (£ = 10). The corresponding graphs show the droplets after the 
initial relaxation (dashed lines, t/C « 60) and in the migration process at x = 36.5 
(solid lines, t/C = 78500, 37500, and 12700 for I, II, and III, respectively). The time 
evolution of the center of mass (x, y) relative to the corner is shown in (b) and (c) for 
the lateral and vertical direction, respectively. 
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Figure 14. The lateral position x of the center of mass as a function of time for 
droplets of size a — 15 (solid line) and a — 10 (dashed line) near the corner of a wedge 
in the minus case with B = 0, C = 2.6667, and 9 eq = 90°. 




Figure 15. (a) Droplets of initial size a = 15 on substrates with 6 eg = 90° (C = 4.2327 
and B = —2.5) starting at I = 10 on the base of steps with heights varying between 
h = 2.5 (top) and h = oo (bottom, corresponding to an isolated wedge) right after the 
initial relaxation process at t = 210 (dashed lines) and during the migration process 
at t = 9100 (solid lines). The horizontal position x and the vertical position y of the 
center of mass are shown in (b) and (c), respectively, as a function of time. 
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side of steps in the minus case and which is more pronounced for higher steps. 

Changing B and C such that 9 eq = 90° remains the same does not change the 
dynamics of the droplets qualitatively. Higher values of B and C result in larger droplet 
velocities (see Fig. 16 concerning the example of droplets near an isolated wedge). Beside 
this change of droplet speed we find that for larger values of C the final position of the 
droplets is closer to the step and the final height y of their center of mass is larger. 
The magnitude of the disjoining pressure and therefore the forces acting on the droplet 
increase with C. In response to these forces the droplets deform more upon increasing 
C. 



5.4- Nanodroplets on edges, wedges, and steps 

In the previous subsections we have discussed the behavior of nanodroplets originally 
positioned at a certain lateral distance from topographic features such as edges, wedges, 
and steps. Their behavior suggests that these surface features provide migration barriers 
for droplets. Even in those situations in which droplets migrate towards the edge or 
wedge, respectively, they stop just before reaching them. This result is also borne out 
in a macroscopic model which takes into account only interface energies: the free energy 
of a droplet positioned right on an edge is larger than that of a droplet of equal volume 
residing on a flat and homogeneous substrate, and the free energy of a drop in the corner 
of a wedge is even lower. As a consequence, we expect that droplets sitting on edges 
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Figure 16. Droplets of initial size a = 15 near isolated wedges (h = oo) in the plus 
case. 6 eq = 90° on both substrates: (I: B = -2.5, C = 4.2327) and (II: B = -4, 
C = 0.9265). (a) The top panel shows the initial droplet shape and the lower panels 
show the droplets just after the initial relaxation at t/C ~ 60 (dashed lines) and in their 
final configuration at t/C = 2500 and t/C = 14500 for substrate I and II, respectively 
(solid lines). The horizontal position x and the vertical position y of the center of the 
mass during the motion are shown in (b) and (c), respectively. 
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to be in an unstable and droplets sitting inside the corner of a wedge to be in a stable 
configuration. Moving, by force, a droplet (with the shape of the liquid- vapor interface 
remaining a part of a circle) in the first case slightly to one side results in an increased 
contact angle on this side, while the contact angle on the other side decreases. However, 
with the leading contact angle being larger than the equilibrium one the corresponding 
contact line will move away from the edge, while the trailing contact line (with the 
corresponding contact angle being smaller than the equilibrium one) moves towards 
the edge. As a consequence, the droplet leaves its position at the edge. In the case 
of a droplet in a wedge, the situation is reversed: moving, by force, the droplet into 
one direction results in a decreased contact angle on this side and an increased contact 
angle at the trailing side, such that the droplet moves back into the corner of the wedge. 
Accordingly one expects that a certain force has to be applied to push a droplet over an 
edge or to pull it out of the corner of a wedge. In the following our detailed numerical 
results indicate that this also holds for nanodroplets and they enable us to quantify 
those external forces. 

Our analyses show that a nanodroplet positioned symmetrically on the tip of an 
edge is unstable on all types of substrates, regardless of whether the droplets migrate 
towards the edge or away from the edge (see Fig. 17; there a suitable, highly symmetric 
initial shape of the liquid-vapor interface has been chosen such that it indeed relaxes 
to the unstable state of a droplet sitting on the tip of the edge). Due to the mirror 
symmetry with respect to the diagonal of the edge, a droplet right at the tip of an 
isolated edge is in mechanical equilibrium but in an unstable one. In the minus case, 
after a small perturbation the droplet flips either up or down but then rests next to 
the step, i.e., in the position which it would assumes upon migrating towards the edge 
(Fig. 17(a)). In the plus expected from the previous results the droplet migrates 

away from the edge after flipping to either side (Fig. 17(b)). At steps of finite height, 
this symmetry is broken by the presence of the wedge. In the minus case, the droplets 
are pushed away from the wedge, i.e., upwards, which is consistent with the dynamics 
of droplets in the vicinity of isolated wedges of the same material. However, being 
attracted to the edge as shown in the previous subsections, they come to rest with the 
trailing contact line pinned to the step edge (Fig. 17(c)). In the plus case the droplets 
move in the opposite direction, i.e., they are attracted by the wedge which they migrate 
to after leaving the edge area (Fig. 17(d)). 

In order to displace a droplet from one side of an edge to the other side (as shown 
in Fig. 18(a)) one has to apply an external force, e.g., a body force such as gravity 
or centrifugal forces, which we incorporate into the hydrodynamic equations via the 
boundary condition (see Eq. (3)). If the applied force is small the droplets assume a new 
but distorted equilibrium position with the leading three-phase contact line still pinned 
at the edge. But there exists a threshold force density gth beyond which the configuration 
described above is unstable and the leading three-phase contact line depins from the 
step edge. As a consequence the droplet flips around the corner and ends up on the 
vertical side of the edge. Since the applied body force has no component parallel to this 
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vertical part of the substrate, the further fate of the droplet is determined by the action 
of the intermolecular forces. In the minus case considered in Fig. 18(a) the droplet is 
attracted to the edge such that the new stable equilibrium configuration is that of a 
droplet residing on the vertical part of the step with the trailing three-phase contact line 
pinned at the step edge. In the plus case (which we have not tested numerically) the 
droplet is repelled from the edge and it is expected to move down the vertical part of 
the edge. As shown in Fig. 18(b) we have determined the body force density g needed 
to push the droplets over the edge for various types of substrates (minus case with 
B — 0) and for droplets of two different sizes. The threshold force density g t h decreases 
both with C (i.e., with 9 eq ) and with the droplet size. Both trends are also expected 
to occur for macroscopic droplets. In the limit 9 eq — > 180° the droplets loose contact 
with the substrate and the free energy of the droplet at the edge equals the free energy 
on a planar substrate. Taking, however, the finite range of molecular interactions into 
account this no longer holds, but the barrier still decreases with increasing 6 eq . Since the 
force density g is a body force, i.e., a force density, the total force per unit ridge length 
G = g Ad (with the ridge cross-sectional area = Q^/ 'L ~ a 2 ) acting on the droplet 
is proportional to the droplet cross-sectional area A^. In Fig. 18(b) we observe that 
the total threshold force Gth needed to push droplets over the edge increases with the 
droplet volume. Apart from the effects of long-ranged intermolecular forces the main 
contribution to the barrier effect of the edge is the increase of the liquid- vapor surface 
area when the droplet is deformed as it passes over the edge. The square root of the 
ratio of the surface tension coefficient 7 and the body force density g defines a capillary 
length below which the surface tension dominates, while it is less important for larger 
drops. (The surface area of three-dimensional droplets increases only quadratically with 
the droplet radius a while the volume increases ~ a 3 .) From dimensional arguments 
gth a 3 ~ 7 a 2 one expects the threshold body force density g t h needed to push droplets 
over an edge to decrease ~ 1/a with the droplet radius, while for liquid ridges the 
total force per unit length G t h should still increase linearly with the droplet radius a. 
Therefore the total threshold force Gth needed for the larger droplet in Fig. 18(b) should 
be about a/288/127 ~ 1.5 times the force needed for the smaller drop. The actual value 
is somewhat smaller and we attribute the difference to the effect of the long-ranged part 
of the intermolecular forces. 

A macroscopic droplet spanning the whole topographic step (i.e., with one contact 
line on the top terrace and one on the base terrace) moves downhill: since the surface 
of a macroscopic droplet is a part of a circle which is cut by the top side of the step 
at a higher level than by the base of the step the contact angle at the upper terrace 
is smaller than the contact angle at the lower terrace. This results in a net driving 
force in downhill direction. The final configuration is a droplet with the upper contact 
line pinned at the step edge. This is also true on the nano scale, as demonstrated in 
Figs. 19(a) and 19(b) for the minus case and in Figs. 19(d) and 19(e) for the plus case. 
The latter indicates that the difference in contact angle at the two contact lines due 
to the different height level at the top side and at the base side of the step provides 
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Figure 17. A droplet positioned symmetrically on an edge (dashed line) is in an 
unstable equilibrium both in the (a) minus case (C = 2.67 and B = 0) and in the (b) 
plus case (C = 4.2327 and B = —2.5). A tiny perturbation at t = pushes the droplet 
cither up or down the step. After that in the minus case the droplet stays next to the 
step (a) whereas in the plus case (b) it moves away from the edge (solid lines, t = 1000 
and t = 1050 in (a) and (b), respectively). In the presence of a wedge, i.e., for a finite 
step height h, the droplet is (c) pushed onto the top side of the step in the minus case 
but (d) towards the corner of the step in the plus case (solid lines, t = 800 and t = 815 
in (c) and (d), respectively). 

a stronger driving force than the lateral action of the disjoining pressure which, in the 
plus case, moves droplets positioned next to the step in uphill direction. 

For all substrates the surface of droplets pinned at the edge becomes convex 
(corresponding to a negative pressure in the droplet) for small 9 eq (i.e., small C, see 
C — 0.1 in Figs. 19(c) and 19(f)). For very large 9 eq the upper contact line depins from 
the edge and moves down towards the wedge (see C = 6 in Figs. 19(c) and 19(f)). The 
result is a droplet sitting in the corner of the wedge area only. The critical value for 
9 eq between both types of configurations depends on the droplet volume and the step 
height: the smaller the droplet (as compared to the step height) the smaller is the value 
of 9 eq at which the upper contact line depins and the larger the volume the smaller is 
the value of 9 eq at which the droplet surface becomes convex. Both phenomena are in 
qualitative agreement with macroscopic considerations which take into account interface 
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Figure 18. (a) A nanodroplet of radius a = 15 pushed over an edge (minus 
case, B = and C = 3) by an external, horizontal body force ge x = 0.00208 e x 
(G = g A c i = 0.71475, direction indicated by the horizontal arrow). Droplet shapes for 
t = (indicated), 100, 650, 1050, 1150, and 6950 are shown (from the upper left to 
the lower right), (b) The minimum total force per unit length Gth = 9th Ad (Ad is 
the droplet cross-sectional area) to push droplets over the edge for two droplet cross- 
sectional areas Ad = 127 (dashed line) and 288 (solid line) corresponding to a ~ 10 
and 15, respectively. The values of {C,B) = (1,0), (2,0), (3,0), and (4,0), correspond 
to 9 eq = 51.3° , 75.5°, 97.2°, and 120°, respectively. The force density (force per unit 
volume) g is measured in units of j/b 2 . 



energies only (see Refs. [25, 68]). 

Droplets sitting in the corner of a wedge are in an energetically rather favorable 
situation as illustrated by the arguments given at the beginning of this subsection. 
However, even in the plus case, for which droplets are attracted by wedges, they stop 
before reaching the wedge and they do not move into the corner of the wedge. In any 
case, there is an energy barrier to overcome in order to move droplets out of wedges, 
as shown in Fig. 20(a). If a small horizontal force is applied to a droplet sitting in the 
corner of a wedge it assumes a new, slightly distorted but stable shape. But there exists 
a threshold force density g t h above which the distorted configuration becomes unstable 
and the droplet moves out of the corner. In the minus case considered in Fig. 20, the 
droplet is repelled from the wedge such that the effect of the intermolecular forces adds 
to the external driving force and the droplet definitively moves out of the corner of the 
wedge. In contrast to the force required to push droplets over an edge, g t h increases 
with C (i.e., with 6 eq ). The total force per unit ridge length G t h = 9th -^d required 
to pull a droplet out of a wedge increases only slightly (i.e., less than linearly) with 
the droplet volume (see Fig. 20(b)), so that accordingly the required force density g th 
decreases significantly with volume. With the same dimensional arguments used above 
for droplets being pushed over edges one would expect the total force needed to pull 
two droplets of different volume out of the corner of a wedge to be proportional to the 
square root of the volume ratio. In particular for small C ~ 2 according to Fig. 20(b) the 
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Figure 19. The final positions and shapes of droplets initially (dashed lines) spanning 
a topographic step of height h = 5 (Ad = 250) and h = 2.5 (Ad = 275) ((a) and (d)), 
h=15 (A d = 150) and h = 10 (A d = 200) ((b) and (e)), and h = 50 (A d = 1667, (c) 
and (f)) in the minus case (B = 0, (a), (b), and (c)) and in the plus case (B = —2.5, (d), 
(e), and (f)), for various values of C . In the minus case C = 0.1, 1, and 4 correspond 
to 6 eq = 15.7°, 51.3°, and 120°, respectively, while for the plus case C = 0.1, 1, 2, and 
6 correspond to 9 eq = 12.5°, 40.2°, 58.2°, and 114.7°, respectively. 
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Figure 20. (a) A droplet with cross-sectional area = 312.5 pulled out of the corner 
of a wedge (B = and C — 3) by an external horizontal body force g e x = 0.00465 
(G = g Ad = 1.453125, direction indicated by the horizontal arrow). Shown are droplet 
shapes for t = (indicated), 400, 1400, 3700, 4700, and 4900 (from left to right), (b) 
The minimum total force per unit ridge length G t h = 9th Ad required to extract the 
droplet from the corner as a function of C for droplets of cross-sectional area Ad = 200 
(corresponding to a ~ 17.5, dashed line) and Ad = 312.5 (corresponding to a ~ 22, 
solid line) in the minus case (B — 0). The values C = 1, 2, 3, and 4 correspond to 
9 eq = 51.3° , 75.5°, 97.2°, and 120°, respectively. The force density g is measured in 
units of 7/6 2 . 



total threshold force is almost independent of the droplet size, rather than to increase 
by a factor a/312.5/200 = 1.25. We attribute this difference to the influence of the 
long-ranged part of the intermolecular forces. 

6. Discussion 

6.1. Force analysis 

Numerical solutions of the Stokes dynamics of nanodroplets in the vicinity of edges, 
wedges, and steps are rather time consuming, even when using advanced numerical 
methods. As shown in Fig. 6(b), the main driving force for the migration of droplets is 
the disjoining pressure induced force density fu as defined in Eq. (6). After the initial 
relaxation process, the shapes of the droplets hardly change during the migration process 
until the droplets either reach the edge (minus case) or the corner of the wedge area 
(plus case). Unfortunately the relaxed shape of the droplet is not available analytically, 
but for droplets on substrates with 9 eq w 90° as mostly considered here the initial shape 
relaxations are rather mild. Accordingly, as demonstrated in the following, the force on 
the droplets can be estimated rather accurately from calculating /n for droplets with a 
shape given by the initial profile in Eq. (5) positioned at the distance t = \x\ — a from 
the edge or from the corner of the wedge. Apparently this estimate becomes invalid for 

Figures 21(a) and 21(b) show the disjoining pressure induced force densities 
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calculated along these lines for droplets of size a = 15 and 5 as a function of the distance 
£ of the right contact line to an edge for the minus and the plus case, respectively. 
Since fu is proportional to C only results for C = 1 are shown. For the plus case 
the force is always negative for both droplet sizes and at all distances from the edge, 
with its strength increasing towards the edge. This means that droplets should move 
away from the edge with a speed which decreases continuously. This is in complete 
agreement with the numerical results presented in the previous section. For the minus 
case and for sufficiently large values of B, the force is positive in accordance with the 
numerical results. However, as shown in Fig. 21(a), for very small values of B, i.e., 
for B < B c ~ — 10 the force in the direct vicinity of the edge becomes negative and 
droplets are expected to move away from the edge. Indeed, as demonstrated in Fig. 8(b) 
(x(t — > oo) for I lies below the corresponding values for II and III), the final distance of 
the droplets from the step edge in the minus case increases with more negative values 
of B. On the other hand, as shown in the following Subsec. 6.2, for large distances from 
the edge, in the minus case even for arbitrarily small B the force is positive so that 
droplets find an equilibrium position with vanishing force at a significant distance from 
the edge. The sign of the disjoining pressure induced force density does not depend on 
the droplet size. However, the equilibrium position changes as a function of droplet size. 

The force calculated for droplets of the same size but in the vicinity of a wedge 
for the minus and the plus case are shown in Figs. 21(c) and 21(d), respectively. For 
the minus case the force is positive for any droplet size and for any B, which means 
that the droplets move away from the wedge. For the plus case the force is negative 
at large distances, but it changes sign close to the wedge at a distance which increases 
with decreasing the size of the droplets and with decreasing the value of B. The latter 
relation is in agreement with the numerical results presented in Fig. 16(b). 

The disjoining pressure induced force density fu presented in Fig. 21 has beeen 
calculated for droplets with a shape given by Eq. (5), i.e., for droplets with equal height 
and half width. However, the substrate parameters used in Fig. 21 do not necessarily 
lead to 6 eq = 90°, and droplets would adopt a very different shape even during the 
migration process. In order to check the influence of the droplet shape on the calculated 
disjoining pressure induced force density fu we also consider droplets which have a width 
w different from their height a (compare with Eq. (5)): 



Figure 22 compares the disjoining pressure induced force density fu on droplets in the 
vicinity of edges and wedges in both the minus and the plus case for different drop 
widths w but for a fixed drop height a = 15. The results indicate that the form of the 
droplets does not change the sign of fu, and in particular in the vicinity of the wedge 
the droplet width has a rather small influence on the force. 

In the vicinity of topographic steps the dependence of fu on the step height h is also 
in good agreement with the results of the full numerical solution of the Stokes dynamics. 




(27) 



CONTENTS 



34 




Figure 21. The DJP induced force density fn (in units of j/b 2 ) on droplets of size 
a = 15 and a = 5 in the vicinity of an edge ((a) and (b)) and a wedge ((c) and (d)) on 
substrates of the minus ((a) and (c)) and the plus ((b) and (d)) type with C = 1 and 
various values of B as indicated in the boxes as a function of the distance I from the 
edge or the corner of the wedge. 



Figure 23 shows fn above and below the step on substrates of the minus and plus type 
for step heights ranging from h = 2.5 to oo (i.e., to isolated edges and wedges). The 
absolute value of the force increases with the step height with the force near isolated 
edges and wedges as the limiting values. This limiting value is almost reached for a step 
height h = 20 (not shown in Fig. 23). 

6.2. Direction of motion far from the step 

Both the force calculations presented in the previous subsection as well as the results 
of the numerical solution of the mesoscopic hydrodynamic equations indicate, that the 
direction of motion of a nanodroplet far enough from the step does not depend on 
whether the droplet is positioned on the top side or on the bottom side of the step. In 
the minus case the droplets move in downhill direction (i.e., in the direction of positive 
x-values) and in the plus case in the opposite direction, independent of the step height 
and of the values of B and C (where the latter has to be positive). In order to understand 
this we further analyze the total force per unit ridge length F n = fn A d on liquid ridges 
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Figure 22. The DJP induced force density /n (in units of 'y/b 2 ) on droplets of height 
a = 15 and widths w = 2 a and u; = 0.5 a in the vicinity of an edge ((a) and (b)) and 
a wedge ((c) and (d)) on substrates of the minus ((a) and (c)) and the plus ((b) and 
(d)) type with C = 1 and various values of B as indicated in the boxes as a function 
of the distance £ from the edge or wedge. 



as defined in Eq. (6) for large droplets far from the step. Asymptotically for large \x\ 
the DJP reduces to its value on a flat substrate so that there the wetting film thickness 
assumes its value yo independent of x up to 0(|x|~ 3 ). The leading order correction to the 
DJP is ±sign(x) 9hC/(16x A ) for the plus and minus case, respectively. Parameterizing 
the shape of the liquid-vapor interface to the left and to the right of the droplet apex 
by xe(y) and x r (y), respectively, from Eq. (6) with dy = —n x ds we obtain: 



Urn 



n(xi{y),y)dy- / U(x r (y),y)dy 
yo Jyo 
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with the droplet apex height y m . In the last but one step we have approximated 



F(x r )-F(xi) 



i x T 
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The force is proportional to the droplet 



cross-sectional area and its sign is determined by the sign of the Hamaker constant (i.e. 
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Figure 23. The DJP induced force density fn (in units of j/b 2 ) on droplets of size 
a = w = 15 on the top side ((a) and (b)) and on the bottom side ((c) and (d)) of steps 
of various heights h with substrates of the minus (B = — 1, (a) and (c)) and the plus 
(B = —2.5, (b) and (d)) type with C = 1 as a function of the distance £ from the step. 



depending on the case; plus or minus) only: in the plus case the force is negative (upper 
sign) and in the minus case it is positive (lower sign). This is in complete agreement 
with the numerical data. However, other than suggested by Eq. (28), the force on a 
droplet does not diverge in the limit h — > oo as this limit has to be taken before taking 
the limit |x| — > oo. 

At large distances from an isolated wedge as well as from an isolated edge, 
the disjoining pressure is to leading order given by the DJP of the corresponding 
homogeneous substrate with ± C/(2 x 3 ) as the leading order correction for the plus and 
the minus case, respectively. As in the case of the step, up to this order the thickness of 
the wetting film is independent of the (large) distances from the edge or wedge. Using 
the same approximations as in the case of the step of finite height, the force on a droplet 
at a distance \x\ from an edge is given by 

Fu = T^- 4 A d , (29) 

with the upper sign corresponding to the plus case and the lower sign to the minus case. 
The sign of the force is the same as in the case of a step and it is also proportional to 
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the droplet volume. However, it decreases less rapidly with the distance from the step. 

For very large, almost macroscopic droplets, the situation is again different from the 
previous two. In the following we follow the line of arguments developed in Refs. [54, 55] 
for droplets in the vicinity of chemical steps. In this limit the droplets are approximately 
symmetric with respect to their apex and the main contribution to the force stems from 
the vicinity of the contact lines. For the wetting film as well as near the apex the x- 
component n x of the surface normal vector is zero and thus in the vicinity of the apex 
the DJP is negligibly small. In most of the examples discussed here the equilibrium 
contact angle 9 eq is about 90° and, as a consequence, the lateral width of the contact 
lines (i.e., the range of x- values within which the drop profile crosses over to the flat one 
of the wetting film) is small and the lateral variation of the DJP within this region is 
negligible. Therefore, after parameterizing the droplet surface in the vicinity of the left 
and right contact line (at x = x — a and x = x + a, respectively) by the corresponding 
function x(y), the total force on a droplet can be approximated by 



roc poo 

/ U(x-a,y)dy- U(x + a,y)dy 

Jyo J m 



= -<$>(x-a,y ) + <$>(x + a,y )xi2ad x <$>(x,yo), (30) 

where $(a;, yo) is the local effective interface potential at the level yo of the wetting 
film (on the top side of the step one has to add h to y ). Extending Eq. (23) to 
inhomogeneous substrates one can define a spatially varying "equilibrium contact angle" 
cos^ e g(x) = 1 + $(x, yo). In this sense, a droplet in the vicinity of a topographic step 
is exposed to an effective chemical wettability gradient which it follows. Expanding 
d x Q(x, y) for large \x\ yields 

. 3 a C h , c ^ f . 

F n ~T—^ + 0(x- 5 ), (31) 

with the upper sign corresponding to the plus case and the lower sign to the minus case. 
The force is equal on both sides of the step and it increases linearly with the step height 
but it decreases rather rapidly with the distance from the step, however, more slowly 
than in the case of nanodroplets. The force increases linearly with the base width 2 a 
rather than with the cross-sectional area Ad- As in the case of nanodroplets the actual 
force on a droplet does not diverge in the limit h — > oo. Using the same approximations 
as in the case of the step of finite height, the force on a droplet at a distance \x\ from 
an isolated edge is given by 

Fn-T^ + C{x~% (32) 

with the upper sign corresponding to the plus case and the lower sign to the minus case. 
In the vicinity of a wedge the situation is more complicated. In order to obtain as in 
Eq. (30) the effective interface potential, the point (x, y) corresponding to the upper 
limit of the integral there has to correspond to a point at infinite distance from the 
substrate. However, taking y to oo for a fixed value x does not change the distance from 
the vertical part of the wedge. At this point it is not clear whether the force integral in 
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Figure 24. The total disjoining pressure induced force per unit ridge length Fn = 
fn Ad at large distance from steps depends on two length ratios h/x and y m /x (the 
cross-sectional area Ad is proportional to y^). This figure summarizes the analytical 
results obtained in Subsec. 6.2. 



Eq. (6) can be approximated by the form given in Eq. (30) because the basic assumption, 
that the disjoining pressure is negligible at the apex, is probably not true. Expanding 
the force as calculated from Eq. (30) for large distances from a step of very large height 
one obtains terms of the order 0(x~ 3 ) competing with terms of order (9(/ix~ 4 ), which 
indicates that in the case of a wedge the approximations Eq. (30) is based on lead to a 
mathematically ill-posed problem. 

In all cases Fu essentially depends on the ratio of the step height h and the distance 
from the step x as well as on the ratio of the apex height y m and x. The asymptotic 
results are summarized in Fig. 24. Fjj varies according to a power law x~^, ( G IN, of the 
distance from the step. For finite sized droplet and steps of finite height (h/x — > and 
Vm/z ~ > 0) we obtain the fastest decay with £ = 5. For almost macroscopic droplets 
{ym/x —> oo) in the vicinity of finite sized steps and for nanodroplets near isolated 
edges and wedges one has £ = 4. For large drops (y m /x — > oo) next to an isolated edge 
we get the weakest decay with £ = 3. In any case, the total force per unit length Fn 
is proportional to the Hamaker constant as observed in the numerical solution of the 
mesoscopic Stokes dynamics as well as in the force analysis presented in Subec. 6.1. 

6.3. Estimates for the velocity 

The driving force / on the droplets is balanced by viscous forces. By applying a simple 
analysis within the lubrication approximation one can show that the rate of energy 
dissipation is proportional to the square of the velocity u = dtx(t) of the droplets 
[1, 69]. The form of this dependence can be expected to hold also for droplets with large 
contact angles on the basis of analyticity and symmetry arguments. By equating this 
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dissipation with the work u f Vt^ done by the driving force one finds 

u = dx/dt ~ f/Q d . (33) 

For droplets far from the step with f(x) given by Eqs. (31) and (32) as a power law one 
has x(t) u ~ t and therefore 

\x(t)\ = (\x \ v + ct) 1 ' v (34) 

with v = 4 for large droplets in the vicinity of edges, v = 5 for large droplets in the 
vicinity of steps of finite height as well as for nanodroplets near isolated edges and 
wedges, and v — 6 for nanodroplets near steps of finite height, c is a constant which 
also depends on whether there is an edge, wedge, or a step and whether the droplet 
is large or small. The functional form given by Eq. (34) with the corresponding value 
of v can be fitted to the positions of nanodroplets as a function of time obtained by 
numerically solving the mesoscopic hydrodynamic equations, e.g., to the data shown in 
Fig. 9 (droplet on the top side of steps, plus case) and Fig. 11 (droplet on the step base, 
minus case). However, the numerically available range of x values is rather small so 
that the fits are consistent with the above values of v but cannot rule out different ones. 
In addition, it is not clear whether the distances considered in the numerical solutions 
of the mesoscopic hydrodynamic equations are large enough to reach the asymptotic 
regime considered here, and whether the droplets should be considered small or large in 
the above sense. 



7. Perspectives 

A major and obvious driving force for studying the dynamics of nanodroplets on 
structured substrates is the rapid development and miniaturization of microfluidic 
devices, in particular of open microfluidic devices [19, 20, 21, 22]. But this is not the only 
research area for which a detailed understanding of the influence of the long-ranged part 
of the intermolecular interactions on fluids in the vicinity of lateral surface structures 
might be important. Another example is the dynamics of nanodroplets at chemical 
surface structures as discussed in Refs. [54, 55]. Moreover, dewetting processes are 
also strongly influenced by surface heterogeneities, both during the initial phase of film 
breakup [16, 70] as well as during hole growth [33]. The latter example is particularly 
interesting in this respect because it reveals an intrinsic nanoscopic length scale which 
has to be understood: the receding contact line is pinned only by steps of a minimum 
height which increases with the size of the liquid molecules [33] — a clear indication that 
details of the intermolecular interactions in the vicinity of the step are relevant. 

All results presented in this article have been obtained for homogeneous straight 
liquid ridges. Apart from the fact that such ridges are unstable with respect to breaking 
up into three-dimensional droplets [71, 72, 23, 24], the question remains to check how 
relevant these results are for actual three-dimensional droplets. In this context we point 
out that the basic driving mechanism for droplets in the vicinity of steps is the difference 
of the disjoining pressure on that side of the droplet which is closer to the step and the 
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side which is further away from the step. In such a situation also three-dimensional 
droplets move. However, the third dimension certainly changes the behavior of droplets 
spanning topographic steps [73, 74]: depending on the droplet volume, the droplet can 
spread along the step into a cigar shaped configuration. The influence of the long-ranged 
part of the intermolecular interactions on this phenomenon has not yet been studied. 

It is worthwhile to point out that, although the dynamics is different, there are 
strong similarities between nanodroplets and solid nanoclusters: their energetics on 
structured surfaces is determined by intermolecular forces as demonstrated in Ref. [75] 
by molecular dynamics simulations of gold clusters on graphite surfaces. Unfortunately, 
in such simulations taking into account the long-ranged component of the intermolecular 
forces increases the numerical cost drastically, such that most of the effects discussed 
here are not accessible by molecular dynamics simulations [76]. 

This leaves the question of experimental tests of our theoretical predictions 
presented here. As detailed in Ref. [56] the forces on the nanodroplets are of the order 
of 10~ 13 N (i.e., about eight orders of magnitude stronger than the gravitational force 
on such a droplet) and the resulting velocities range between O.lmm/s and 0.1/xm/s 
for viscosities between 0.1 Pas and 100 Pas. While topographic surface structures of 
almost any type can be produced with modern lithographic techniques, positioning 
nanodroplets with nanometer accuracy next to a step remains a tough challenge. Most 
promising are techniques based on using atomic force microscopes as pens [77, 78], 
but there are no experiments available yet. Experimentally it is much easier to grow 
droplets from an aerosol or a vapor phase rather than to deposit them at a specific 
location. Experimentally it has been shown, that water nanodroplets preferentially 
condense onto terrace steps of vicinal surfaces [79]. However, these data do not allow 
one to determine whether the droplets reside on the top terrace, on the bottom terrace, or 
whether they span the step. This example shows, that condensation on (and evaporation 
from) nano-structured substrates is a challenging problem of its own. 
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Appendix A. Numerical algorithm 

In order to study the effect of the intermolecular forces on nanodroplets near a 
topographic step we solve Eqs. (l)-(3) numerically using a standard and accurate 
biharmonic boundary integral method (BBIM) [80, 81, 82, 56, 54, 55]. To this end we 
introduce the stream function i/j(x, y) so that dip/dy = u x and dip/dx = —u y as well as 
the vorticity u(x, y) = du x /dy—du y /dx which allows us to reformulate the dimensionless 
versions of Eqs. (1) and (2) in terms the following harmonic and biharmonic equations 




i /+ I 
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Figure Al. The dynamics of a droplet is investigated which initially is positioned 
either on the top side or the bottom side of a step at a distance £ from the step. 
The boundary of the system T decomposes into three different groups: liquid-liquid 
(11), liquid-solid (Is), and liquid- vapor (Iv) interfaces. The discretized node points 
considered in the numerical investigation are indicated; n and t represents the normal 
and the tangential unit vectors on T, respectively. 



[80, 81, 82, 56, 54, 55]: 

V 2 w = (A.l) 

and 

VV = • (A.2) 

The standard BBIM relies on mapping the equations for uj and ip onto the boundary 
r(s) = (x(s),y(s)) of the fluid, parameterized in terms of its contour length parameter 
s. This results in an integral equation for u, ip, and their derivatives u n = n • Vc<j 
and ip n = n • V^>, with the surface normal vector pointing outwards of the liquid. 
By dividing the boundary of the system into a series of elements (see Fig. Al) one 
obtains a coupled system of algebraic equations which can be solved numerically. With 
the tangential velocity ut = ip n and the normal velocity u n = —ip s (with the index s 
indicating the derivative in the direction tangential to the boundary) the position of 
the liquid boundary after a time step can be calculated via the explicit Euler scheme 
[80, 81, 82, 56, 54, 55]: 

r(t + At) = r(t) +u(t)At (A.3) 

In order to solve these equations the boundary conditions of the system must be 
expressed in terms of uo and ip. Depending on the phases in contact with each other, 
three different types of boundary interfaces can be identified (see Fig. Al): liquid-solid 
interfaces T/ s , liquid-liquid interfaces (those boundaries Vm and r« 2 which are located 



CONTENTS 



42 



at the end sides of the system), and liquid- vapor interfaces Ti v . For Ti s we impose the 
no-slip condition (u = 0) which corresponds to ip = and i/) n = 0. For Tm and Tu 2 
we apply a no-flux condition which corresponds to having a vertical symmetry plane 
there. Such a system corresponds to a periodic repetition of the system attached to its 
mirror image. Correspondingly the slope of the liquid-vapor interface at the side ends 
of the system is zero. These conditions can be implemented by setting ip = and u = 
there. The tangential and the normal component of the boundary condition (3) along 
the liquid-vapor interface Ti v in terms of the stream function and the vorticity read 
(lower indices s indicate derivatives with respect to the contour length parameter s) 

u) = 24) ss + 2n4) n (A.4) 

and 

u n = -2 ip nss + 2 k iJ ss + 2 k s tp s + Ks + n ' + - Xs , (A. 5) 



C 



respectively, with the local curvature 



_ y ss x s x ss y s 

W + V*)*» " ( } 

In order to increase the efficiency of the numerical calculations we employ an 

adaptive time stepping: for any numerical step, the time step is selected such that 

the displacement of any node does not exceed 5 percents of the length of the elements 

connected to that node; 5 can be changed during the numerical calculations. The 

starting value for 5 and the rate of its increase depends on the actual situation but 

typically we have started with 5 = 0.01 and then gradually increased it to 0.1 or even 

more. In order to avoid numerical instabilities, the position of the end points of the 

boundary elements are smoothed after a specified number of steps by fitting a spline 

through the points on the liquid-vapor interface, followed by selecting new and equally 

spaced points on the spline. 
[1] dc Gcnncs P G 1985 Rev. Mod. Phys. 57 827-860 

[2] Dietrich S in Phase Transitions and Critical Phenomena, ed. Domb C and Lebowitz J L Vol. 12; 

Academic Press, London 1988; chapter 1, pp. 1-218 
[3] Dietrich S In ed. Caccamo C, Hansen J P and Stell G, New approaches to old and new problems 

in liquid state theory Vol. C529 of NATO-ASI Series Cpp. 197-244 Dordrecht 1999. Kluwer 
[4] Swain P S and Lipowsky R 1998 Langmuir 14 6772-6780 
[5] Richard D and Quere D 1999 Europhys. Lett. 48 286-291 
[6] Lafuma A and Qucrc D 2003 Nature Materials 2 457-460 

[7] Krupenkin T N. Taylor J A, Schneider T M and Yang S 2004 Langmuir 20 3824-3827 
[8] Quere D 2005 Rep. Prog. Phys. 68 2495-2532 

[9] Yang C, Tartaglino U and Persson B N J 2006 Phys. Rev. Lett. 97 116103 
[10] Sbragaglia M, Peters A M. Pirat C, Borkent B M, Lammcrtink R G H, Wessling M and Lohse D 

2007 Phys. Rev. Lett. 99 156001 
[11] Quere D 2008 Ann. Rev. Mater. Sci. 38 71-99 
[12] Bico J, Tordcux C and Qucrc D 2001 Europhys. Lett. 55 214-220 

[13] McHale G, Shirtcliffc N J, Aqil S, Perry C C and Newton M I 2004 Phys. Rev. Lett. 93 036102 
[14] Extrand C W, Moon S I. Hall P and Schmidt D 2007 Langmuir 23 8882-8890 
[15] Rockford L, Liu Y, Mansky P, Russell T P, Yoon M and Mochric S G J 1999 Phys. Rev. Lett. 82 
2602-2605 



CONTENTS 



13 



[16] Kargupta K and Sharma A 2002 Langmuir 18 1893-1903 

[17] Mitchell P 2001 Nature Biotech. 19 717-721 

[18] Thorsen T, Maerkl S J and Quake S R 2002 Science 298 580-584 

[19] Zhao B, Moore J S and Becbc D J 2001 Science 291 1023-1026 

[20] Lam P, Wynne K J and Wnck G E 2002 Langmuir 18 948-951 

[21] Zhao B, Moore J S and Becbc D J 2002 Anal. Ghem. 74 4259-4268 

[22] Zhao B, Moore J S and Becbc D J 2003 Langmuir 19 1873-1879 

[23] Koplik J, Lo T S, Rauscher M and Dietrich S 2006 Phys. Fluids 18 032104 

[24] Mechkov S, Rauscher M and Dietrich S 2008 Phys. Rev. E 77 061605 

[25] Sccmann R, Brinkmann M, Kramer E J, Lange F F and Lipowsky R 2005 PNAS 102 1848-1852 

[26] Robbins M O, Andclman D and Joanny J F 1991 Phys. Rev. A 43 4344-4354 

[27] Nctz R R and Andelman D 1997 Phys. Rev. E 55 687-700 

[28] Rejmer K, Dietrich S and Napiorkowski M 1999 Phys. Rev. E 60 4027-4042 

[29] Rascon C and Parry A O 2000 Nature 407 986-989 

[30] Bruschi L, Carlin A and Mistura G 2002 Phys. Rev. Lett. 89 166101 

[31] Klicr J, Lcidcrcr P and Wyatt A F G 2005 Phys. Rev. B 72 245410 

[32] Gang O, Alvinc K J, Fukuto M, Pershan P S, Black C T and Ocko B M 2005 Phys. Rev. Lett. 95 
217801 

[33] Ondarcuhu T and Picdnoir A 2005 Nano Letters 5 1744-1750 

[34] Tasinkevych M and Dietrich S 2006 Phys. Rev. Lett. 97 106102 

[35] Tasinkevych M and Dietrich S 2007 Eur. Phys. J. E 23 117-128 

[36] Gramlich C M, Mazouchi A and Homsy G M 2004 Phys. Fluids 16 1660-1667 

[37] Dupuis A and Yeomans J M In ed. Bubak M, van Albada G D, Sloot P M A and Dongarra J J, 

Computational Science - ICCS 2004 Vol. 3039 of Lecture Notes in Computer Science pp. 556-563 

Berlin 2004. Springer 

[38] Dupuis A and Yeomans J M 2005 Langmuir 21 2624-2629 

[39] Davis J M and Troian S M 2005 Phys. Fluids 17 072103 

[40] Gaskcll P H, Jimack P K, Sellier M and Thompson H M 2006 Phys. Fluids 18 013601 

[41] Saprykin S, Koopmans R J and Kalliadasis S 2007 J. Fluid Mech. 578 271-293 

[42] Tseluiko D, Blyth M G, Papageorgiou D T and Vanden-Broeck J M 2008 J. Fluid Mech. 597 
449-475 

[43] Oron A 2000 Phys. Fluids 12 1633-1645 

[44] Eijkel J C T and van den Berg A 2005 Microfluid Nanofluid 1 249-267 

[45] Mukhopadhyay R 2006 Anal. Chem. 78 7379-7382 

[46] Rauscher M and Dietrich S 2008 Ann. Rev. Mater. Sci. 38 143-172 

[47] Fetzer R, Rauscher M, Munch A, Wagner B A and Jacobs K 2006 Europhys. Lett. 75 638-644 

[48] Fetzer R, Rauscher M, Seemann R, Jacobs K and Mcckc K 2007 Phys. Rev. Lett. 99 114503 

[49] Brusch L, Kiihne H, Thick U and Bar M 2002 Phys. Rev. E 66 011602 

[50] Bielarz C and Kalliadasis S 2003 Phys. Fluids 15 2512-2524 

[51] Thiele U, Brusch L, Bestehorn M and Bar M 2003 Eur. Phys. J. E 11 255-271 

[52] Gaskcll P H, Jimack P K, Sellier M and Thompson H M 2004 Int. J. Numer. Meth. Fluids 45 
1161-1186 

[53] Yochclis A, Knobloch E and Pismcn L M 2007 Eur. Phys. J. E 22 41-49 

[54] Moosavi A, Rauscher M and Dietrich S 2008 Langmuir 24 734-742 

[55] Moosavi A, Rauscher M and Dietrich S 2008 J. Chem. Phys. 129 044706 

[56] Moosavi A, Rauscher M and Dietrich S 2006 Phys. Rev. Lett. 97 236101 

[57] Bauer C and Dietrich S 1999 Phys. Rev. E 60 6919-6941 

[58] Lauga E, Brenner M P and Stone H A in Springer Handbook of Experimental Fluid Mechanics, 
ed. Tropea C, Yarin A L and Foss J F; Springer, Berlin 2007; chapter Part C, pp. 1219-1240 

[59] Griin G, Mcckc K R and Rauscher M 2006 J. Stat. Phys. 122 1261-1291 

[60] Chaudhury M K and Whitesides G M 1992 Science 256 1539-1541 



CONTENTS 



44 



[61] Subramanian R S, Moumen N and McLaughlin J B 2005 Langmuir 21 11844-11849 

[62] Pismcn L M and Thick U 2006 Phys. Fluids 18 042104 

[63] Ondarcuhu T and Raphael E 1992 C. R. Acad. Sci. Paris, Serie II 314 453-456 

[64] Dietrich S and Napiorkowski M 1991 Physica A 177 437-442 

[65] Dietrich S and Napiorkowski M 1991 Phys. Rev. A 43 1861-1885 

[66] Napiorkowski M, Koch W and Dietrich S 1992 Phys. Rev. A 45 5760-5770 

[67] Raphael E 1988 C. R. Acad. Sci. Paris, Seme II 306 751 

[68] Kharc K, Hcrminghaus S, Barct J C, Law B M, Brinkmann M and Seemann R 2007 Langmuir 23 
12997-13006 

[69] Scrvantic J and Miiller M 2008 J. Chem. Phys. 128 014709 

[70] Mukherjee R, Pangule R C, Sharma A and Banerjee I 2007 J. Chem. Phys. 127 064703 

[71] Davis S H 1980 J. Fluid Mech. 98 225 

[72] Brinkmann M, Kierfeld J and Lipowsky R 2005 J. Phys.: Condens. Matter 17 2349-2364 

[73] Brinkmann M and Blosscy R 2004 Eur. Phys. J. E 14 79-89 

[74] Herminghaus S, Brinkmann M and Seemann R 2008 Ann. Rev. Mater. Sci. 38 101-121 

[75] Yoon B, Luedtke W D, Gao J and Landman U 2003 J. Phys. Chem. B 107 5882-5891 

[76] De Coninck J and Blake T D 2008 Ann. Rev. Mater. Sci. 38 1-22 

[77] Yang M, Sheehan P E, King W P and Whitman L J 2006 J. Am. Chem. Soc. 128 6774-6775 

[78] Fang A, Dujardin E and Ondarcuhu T 2006 Nano Letters 6 2368-2374 

[79] Hu J, Carpick R W, Salmeron M and Xiao X D 1996 J. Vac. Sci. & Tech. B 14 1341-1343 

[80] Kelmanson M A 1983 J. Eng. Math. 17 329-343 

[81] Betelii S, Diez J, Thomas L, Gratton R and Marino B 1997 Int. J. Numer. Meth. Fluids 25 1-19 

[82] Mazouchi A, Gramlich C M and Homsy G M 2004 Phys. Fluids 16 1647-1659 



